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A  FREE  ELECTRON  LASER  DRIVEN  BY  A 
LONG  PULSE  INDUCTION  LINAC 


I.  Introduction 

In  the  last  five  years,  free  electron  laser  (FEL)  research  has  expanded 

rapidly.  Great  progress  has  been  made  both  in  the  theoretical  understanding 

of  the  mechanism  and  in  the  experimental  demonstration  of  the  validity  of  many 

features  predicted  by  theory.  The  FEL  has  attracted  this  interest  because  of 

its  promise  as  a  tunable,  powerful,  and  efficient  source  of  radiation  over  a 

very  large  frequency  range.  The  strong  dependence  of  output  frequency  on 

2 

electron  energy  (u  <*  Y  )  makes  scaling  attractive.  Furthermore,  the 
availability  of  very  high  power  electron  beams  (>>1  GW)  makes  possible  the 
generation  of  high  power  radiation  in  regions  of  the  spectrum  where  other 
practical  high  power  sources  do  not  exist. 

In  this  paper  we  will  present  an  overview  of  FEL  research  and  then 
describe  in  some  detail  an  ongoing  experiment  designed  to  produce  mm-wave 
output  from  a  long  pulse  induction  linac  (t  *  2  ysec) .  in  Section  II  we 
will  describe  the  FEL  mechanism,  and  a  brief  history  of  FEL  experimental 
research  will  be  presented  in  Section  III.  Section  IV  will  describe  the 
various  FEL  operating  regimes  and  discuss  the  beam  requirements  for  operation 
in  the  Raman  regime.  Accelerator  development  and  diode  characteristics  will 
be  discussed  in  Section  V  and  the  long  pulse  induction  linac  described.  Our 
FEL  experiment  will  be  described  in  Section  VI  along  with  a  comparison  of 
various  wigglers.  The  experimental  results  will  be  presented  in  Section  VII 
and  future  directions  for  the  experiment  and  conclusions  in  Section  VIII. 
Manuscript  approved  December  28, 1982. 
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II.  FEL  Mechanism 


To  understand  the  FEL,  let  us  consider  a  very  low  current  electron  beam  so 
that  space  charge  forces  are  not  important.  We  assume  the  initial  velocity  of 
the  electron  beam  is  in  the  axial  direction  and  we  have  a  simple  linear 
wiggler  in  the  y  direction. 

Let 


and 


A 

B  -  B  sin(k  z)  y  . 
w  w  w  ■* 


Then 


v  «  v  cos  (k  z)  x  . 
w  w  w 


(2) 

(3) 


Here,  vb  is  the  beam  velocity,  Bw  the  wiggler  magnetic  field,  and  vw  is  the 
wiggle  velocity  that  results  from  the  Lorentz  force  acting  on  the  particle  as 
it  passes  through  the  wiggler  field,  we  now  assume  the  presence  of  a  linearly 
polarized  radiation  field 

A  A 

E  +  B  -  E  cos  (kz  -  wt)  x  +  B  sin  (kz  -  wt)  y  . 


This  radiation  field  will  exist  as  part  of  the  noise  spectrum  in  the  case  of 
the  oscillator,  or  will  be  supplied  externally  in  the  case  of  an  amplifier. 
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As  a  result  of  the  wiggler  and  the  radiation  field,  a  ponderraotive  force 


in  the  z  direction  develops 


F  «^(v  xB  +  vxB  )  z 
—  c  w  —  —  w 


sin[(k  +  k  )  z  -  ajt]  z 
-  w 


(5) 


This  force  arises  from  the  interaction  between  vw  and  the  magnetic  cong>onent 
of  the  radiation  field  and  from  the  perturbed  velocity  v  (which  results 
from  v.  x  B)  interacting  with  the  wiggler  field. 

P  “ 

This  ponderraotive  force  drives  a  current,  6 Jz .  From  the  continuity 

equation  q  3<Sn/3t  =  7  •  5j  ,  we  obtain  a  density  modulation  of  the  form 

z 

Sn  *  cos  ( (k  +  k^)  z  -  rat]  .  (6) 

This  density  modulation  of  the  electron  beaus,  driven  by  the  pondermotive 
force,  results  in  a  current  which  can  cause  radiation  in  phase  with  the 
radiation  field: 


6j  »  q  fin  v^  «  cos  (kz  -  o»t)  x  .  (7) 

Therefore  the  pondermotive  force  causes  the  beam  to  bunch  and  the 
electrons  to  radiate  coherently  and  in  phase  with  the  existing  radiation. 
Thus  the  original  radiation  grows,  which  increases  5n,  increasing  $J  and 
so  on  until  saturation  occurs. 


If  An[  (k  +  k  )  ,  u>]  is  at  a  frequency  and  wave  number  that  is  also  a  beam 
mode,  then  we  have  a  collective  free  electron  laser.  Then  to  determine  the 
wavelength  scaling,  we  assume  the  phase  velocity  of  the  pondermotive  wave  is 
near  the  beam  velocity 


ut 

V  *  -  1  "  1  •  v 

ph  k+k  b  * 
w 


(8) 


Then  using  w  ■  ck,  we  obtain 


X 


(9) 


This  is  the  wavelength  scaling  law  that  is  such  an  important  feature  of  the 
FEL.  If  the  beam  current  is  high,  the  pondermotive  force  will  drive  a 
collective  space  charge  wave.  This  enhances  the  interaction  and  increases  the 
intrinsic  efficiency.  However  the  higher  current  introduces  an  additional 
complication,  that  of  the  beam  rapidly  expanding  due  to  its  own  space 
charge.  The  siiqplest  configuration  to  confine  the  beam  is  an  axial  magnetic 
field.  This  introduces  the  cyclotron  modes  on  the  beam,  which  can  result  in 


cosseting  interactions 


'  ’3>*r'^VW‘Cir,J 


III.  fel  Experimental  Research 


As  is  typical  in  a  rapidly  developing,  complex  field  such  as  free  electron 
lasers,  experiments  have  lagged  behind  theory  to  a  large  degree.  In  the  last 
2-3  years,  however,  scores  of  experiments  have  been  proposed  and  undertaken  to 
explore  the  many  aspects  of  FEL  operation  predicted  by  theory.  Although  only 
preliminary  results  have  been  reported  from  most  of  these  experiments, 
indications  are  that  the  FEL  mechanism  is  a  viable,  tunable  source  of 
radiation  from  millimeter  waves  to  the  infrared.  Experiments  now  underway  are 
attempting  to  extend  the  range  of  operation  into  the  visible  and  even  beyond. 


Free  electron  laser  experiments  have  evolved  along  two  distinct  paths 
depending  on  the  type  of  electron  accelerator  used.  High  current,  relatively 
low  particle  energy  electron  beams  typically  from  Marx  type  generators  have 
been  used  for  FEL  experiments  in  the  collective  regime.  More  conventional 
accelerators  such  as  RF  linacs,  which  produce  much  higher  energy  electrons  but 
at  much  lower  current  levels,  have  proved  useful  for  FEL  research  in  the 
Compton  or  single  particle  regime,  Madey  and  coworkers1  at  Stanford  have 
pioneered  the  Compton  regime  research  using  a  superconducting  linac  together 
with  a  helical  magnetic  wiggler  field  (Elias  et  al, ,  1976;  Deacon  et  al. , 
1977).  They  initially  performed  amplifier  experiments  at  X  ■  10.6  ym  and 
later  demonstrated  laser  oscillator  operation  at  X  *  3,4  ym  using  a  43  MeV 
beam  having  a  peak  micropulse  current  of  2.6  A.  Hie  efficiency  (laser 
energy/beam  energy)  of  the  oscillator  was  less  than  10“*,  but  it  was  suggested 
that  a  much  higher  efficiency  could  be  achieved  by  providing  for  multiple 
passes  of  the  electron  beam  through  the  FEL  interaction  region.  However,  the 
increasing  beam  energy  spread  with  each  pass  becomes  a  limiting  factor. 
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More  recently ,  the  Stanford  group  has  collaborated  with  LURE  on  an 
experiment  using  the  ACO  storage  ring  in  Qrsay,  France  (Bazin  et  al. ,  1982). 

In  this  experiment  a  superconducting  undulator  (B^  *  4  kG,  L  ~  1  m)  was 
placed  on  a  straight  section  of  the  storage  ring/  which  was  typically  operated 
at  150  MeV  with  a  bunch  current  of  ^  10  mA.  An  argon 
laser  (X  »  4880  A  or  5145  A)  having  a  power  density  of  ~  1.6  kW/cm2  was 

•  4 

amplified  in  the  interaction  region,  lhe  maximum  gain  reported  was  **•  4  x  10 
per  pass,  hence  oscillator  experiments  would  be  extremely  difficult  with  the 
original  apparatus.  Consequently ,  the  group  is  attempting  to  increase  the 
gain  by  substituting  an  optical  klystron  for  the  undulator. 


Somewhat  similar  experiments  are  underway  at  Frascati  using  the  Adone 
storage  ring  (Barbini  and  Vignola,  1982)  and  at  Brookhaven  using  the  VUV 
storage  ring  (Luccio,  1982) .  Another  storage  ring  FEL  experiment  has  been 
proposed  for  BESSY  in  Berlin  (Gaupp,  1982). 


Another  set  of  similar  experiments  eiqploying  rf  linacs  in  conjunction  with 
C02  lasers  and  planar,  permanent  magnet  wigglers  have  recently  been 
undertaken.  These  experiments  were  initiated  at  Los  Alamos,  TRW,  and  Math 
Sciences  Northwest.  Hie  permanent  magnet  wigglers  used  in  these  experiments 
are  convenient  for  studying  variable  wiggler  efficiency  enhancement  schemes  by 
changing  either  the  wiggler  period  or  amplitude  with  axial  distance.  In 
initial  experiments,  the  groups  have  concentrated  on  measuring  the  energy  loss 
of  the  electrons  passing  through  the  interaction  region,  because  the  gains  or 
efficiencies  are  so  low  that  it  is  very  difficult  to  accurately  measure  the 
amplified  output  signal  in  the  presence  of  the  large  input  laser  signal  (20  - 
l000  MW) . 
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Based  on  measurements  of  electron  energy  loss,  the  TRW  group  has 
calculated  an  efficiency  of  0.07%  and  a  gain  of  2.7%  using  a  wiggler  with  an 
amplitude  taper  of  2.25%  (Boehmer  et.  al. ,  1982).  This  efficiency  is  an  order 
of  magnitude  higher  than  the  theoretical  value  with  no  taper.  By  using  a  more 
powerful  laser  than  the  TRW  group  (1  GW  vs  2  MW)  the  Los  Alamos  group  can 
achieve  high  efficiency  (~  2%)  but  lower  gain  (Warren  et.  al. ,  1982).  (Note 
that  the  input  laser  power  is  not  included  in  the  efficency  calculation) .  The 
Los  Alamos  experiment  is  now  being  converted  to  an  oscillator  mode. 

Additional  low-current-regime  experiments  are  being  performed  using 
microtron  accelerators  at  Bell  Labs  (Shaw  and  Batel,  1982)  and  at  Frascati 
(Bizzarri  et  al. ,  1982).  These  accelerators  operate  at  up  to  ~20  MeV  with 
peak  current  of  ~  5  A,  and  the  Frascati  microtron  is  being  upgraded  to  >  30 
MeV.  Bell  Labs  is  using  a  10  m  long  helical  wlggler  with  a  20  cm  period  to 
generate  100  -  400  Jim  radiation  by  tuning  the  beam  energy  from  10-20  MeV. 
Frascati  uses  a  2.25  m  long  permanent  magnet  wiggler  with  a  period  of  5  cm  for 
output  at  -  15-20  wm. 

A  different,  but  quite  interesting  approach  is  being  pursued  at  the 
university  of  California  -  Santa  Barbara  (Elias  and  Ramian,  1982).  The  UCSB 
group  is  using  a  recirculating  electrostatic  accelerator  to  achieve  high 
average  power,  good  beam  quality,  and  high  overall  FEL  efficiency.  If  all  the 
beam  could  be  recovered,  the  device  could  essentially  operate  dc.  The  group 
hopes  to  use  the  2A,  3  MeV  electron  beam  to  generate  FIR  to  submillimeter 
radiation  in  a  single  stage  FEL,  and  then  to  perform  a  two  stage  experiment  to 
generate  visible  to  IR  output. 
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High  current,  relatively  low  particle  energy  FEL  experiments  have  achieved 
impressively  high  gains  and  radiation  power  levels  at  much  lower  frequencies 
than  the  high  particle  energy  experiments*  However,  these  experiments  are 
characterized  by  a  different  set  of  problems.  High  current  beams  are 
typically  produced  by  Marx-type  generators  and  in  general  have  higher  energy 
spreads  than  low  current  beams.  Self  field  effects  become  important  in  the 
propagation  of  these  beams,  and  the  axial  magnetic  field  requried  to  confine 
the  beam  can  result  in  large  amounts  of  cyclotron  type  emission. 

The  first  operation  of  a  device  employing  what  we  now  call  the  FEL 
mechanism  was  reported  over  20  years  ago  by  Phillips  (1960).  He  called  the 
device  a  Ubitron  and  built  several  microwave  tube-type  versions  which  behaved 
impressively  in  S-band  (~  2-3  GHz).  The  wiggler  typically  consisted  of 
alternately  wound  co-axial  coils  and  the  interaction  length  was  ~  50  cm.  Peak 
microwave  output  power  in  excess  of  1  MW  was  achieved  using  a  beam  of 

150  kV  and  <  100  A,  corresponding  to  an  efficiency  of  10%.  Scaling  this 

device  into  the  ram  wavelength  range  was  restricted  at  the  time  because  of  the 

limited  voltage  of  available  electron  guns.  However,  the  introduction  of  high 

current,  higher  voltage  Marx  type  generators  made  possible  a  renewed  effort  in 
the  generation  of  high  power  mm  waves  using  this  mechanism. 

High  current  (>  1  kA)  FEL  experiments  were  originally  referred  to  as 
stimulated  scattering,  and  their  origins  can  be  traced  to  early  experiments  at 
Cornell  University  and  the  Naval  Research  Lab  in  which  large  radiated  powers 
were  observed  when  an  intense  beam  was  modulated  by  a  slow  wave  structure  or  a 


rippled  magnetic  field  (Nation,  1970;  Friedman  and  Herndon,  1972).  These  and 


similar  experiments  were  interpreted  using  one  of  two  different  theories,  the 

cyclotron  maser  instability  (Sprangle  and  Manheimer,  1975)  and  stimulate 

Raman  scatter ing{ Sprangle  et.  al,  1975),  which  predicted  large  growth  rates 

and  high  efficiencies.  The  cyclotron  maser  instability  is  a  purely 

relativistic  effect  associated  with  relativistically  gyrating  electrons.  This 

mechanism  has  been  successfully  used  in  the  development  of  gyrotrons, 

relatively  contact  tube-type  devices  which  operate  at  the  fundamental  or 

harmonics  of  the  cyclotron  frequency  (See,  e.g.,  Granatstein,  ed. ,  1981). 

Stimulated  scattering  experiments,  on  the  other  hand,  have  continued  to  employ 

Marx  generators  at  higher  voltage  and  current  levels  in  order  to  take 

2 

advantage  of  the  y  frequency  scaling  and  the  growth  rate  scaling  with 
current.  After  the  Stanford  group  popularized  the  term  "free  electron  laser" 
and  the  equivalence  of  the  FEL  and  stimulated  scattering  mechanisms  was  shown 
(Kroll  and  McMullin,  1978;  Sprangle  et.  al. ,  1979),  the  high  current 
experiments  also  began  to  be  referred  to  as  FEL’s. 

The  first  high  current,  relativistic  beam  experiment  to  be  interpreted  as 
stimulated  scattering  was  performed  at  the  Naval  Research  Laboratory 
(Granatstein  et  al. ,  1974).  In  experiments  designed  to  produce  microwaves  at 
X  *  2  cm  via  the  cyclotron  maser  interaction,  strong  submillimeter  radiation 
was  also  observed.  A  quantitatively  consistent  explanation  was  that  some  of 
the  2  cm  radiation  reflected  off  the  output  window  and  subsequently  interacted 
with  a  sufficiently  cold  part  of  the  1.5  MeV  beam  to  produce  the  high 
frequency  scattered  radiation.  In  a  subsequent  experiment  designed  to 
optimize  this  effect,  400  ym  radiation  was  generated  at  a  power  of  ~  1  MW 
(Granatstein  et  al.,  1977). 
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Hie  next  important  step  in  the  development  of  high  current  PEL'S  came  in 
an  experiment  at  Columbia  University  in  which  the  high  power  electromagnetic 
pimp  was  replaced  by  a  periodic  magnetic  wiggler,  or  quasi-wave  (Mross  et  al. , 
1976;  Efthimion  and  Schlesinger,  1977).  A  750  keV,  5  kA  beam  was  used  to 
generate  >  10-100  kW  of  mm  or  cm  radiation.  Subsequent  experiments  at 
Columbia  increased  the  millimeter  wave  power  level  to  more  than  1  MW  (Marshall 
et  al.,  1977a).  This  group  also  performed  mode  structure,  spectral,  and 
growth  rate  measurements  which  indicated  that  the  mechanism  responsible  for 
the  observed  output  was  weak  pump  Raman  scattering  (Marshall  et  al. ,  1977b; 
Gilgenbach  et  al.  ,  1979). 

These  early  e)q>eriments  at  NRL  and  Columbia  led  to  a  collaborative  effort 
on  NRL's  VEBA  generator  in  which  a  quasioptical  cavity  was  used  to  provide 
feedback  (McDermott  et  al. ,  1978).  Approximately  1  MW  of  400  um  radiation  was 
coupled  out  of  the  cavity,  and  the  output  wavelength  agreed  with  theoretical 
predictions  for  the  y  ■  3.4  beam.  An  important  result  of  this  experiment  was 
the  demonstration  of  line  narrowing  when  feedback  was  introduced.  The  line 
width  decreased  from  AX/X  >  10%  to  AX/X  *  2%. 

The  efficiency  (<  0.03%)  of  these  high  current  experiments  was  strongly 
limited  by  the  poor  beam  quality.  Also,  researchers  at  MIT  have  shown  that 
the  output  from  such  experiments  can  easily  be  dominated  by  cyclotron  emission 
(Shefer  and  Bekefi,  1982).  Consequently,  the  NRL  VEBA  group  (Parker  et  al. , 
1982)  designed  an  improved  diode  in  which  90%  of  the  diode  current  was  removed 
from  the  beam  with  a  collimator.  The  collimated  1.5  kA,  1.4  MeV  beam  had  an 
instantaneous  axial  velocity  spread  of  ~  0.1%.  The  beam  was  passed 
through  a  helical  wiggler  and  generated  ~  35  MW  of  radiation  at 


X  *  4  mm.  Itiis  power  was  ~  2.5%  of  the  propagated  beam  power.  A  strong 
frequency  dependence  of  the  output  on  the  axial  guide  field  was  observed,  with 
the  highest  output  occuring  at  a  field  slightly  above  the  cyclotron 
resonance. 

Another  recent  collective  regime  experiment  at  Columbia  University  has 

been  reported  in  which  the  electrons  were  given  an  initial  transverse  velocity 

before  their  injection  into  the  wiggler  (Grossman  et.  al. ,  1983) .  Output  was 

observed  at  a  frequency  equal  to  the  sum  of  the  doppler  shifted  cyclotron  and 

2 

usual  FEL  frequencies;  i.e.,  w  *  2y  (fl  /y  +  k  v  ) .  In  this  experiment, 

o  z  o  w  z 

<  1  MW  of  power  at  X  *  1.5  mm  was  generated.  lhis  mechanism  was 
independently  proposed  by  the  MIT  group,  which  called  it  a  Lowbitron  (McMullin 
and  Bekefi,  1981),  and  experiments  to  study  the  effect  are  also  underway 
there. 

Finally,  an  important  trend  that  is  occurring  in  FEL  research  is  an 
attempt  to  bridge  the  gap  between  the  two  current  regimes.  High  current 
experiments  are  being  scaled  to  higher  voltages  in  an  attempt  to  produce 
higher  output  frequencies  and/or  higher  powers.  In  an  experiment  being 
undertaken  on  the  ETA  induction  linac  (V  ~  4.5  MV)  at  Livermore,  a  high 
current  beam  (  ~  1  kA)  is  being  used  in  an  attempt  to  efficiently  generate 
very  high  FEL  output  powers  at  f  ~  100  GHz  (Prosnitz  and  Sessler,  1982).  At 
the  Naval  Research  Laboratory,  a  program  is  underway  to  develop  compact  high 
current,  and  high  voltage  accelerators,  which  could  be  used  as  FEL  drivers 
(Raberson  et  al.  ,  1982).  Die  long  pulse  duration  experiment' reported  here  is 
part  of  that  program  and  is  scalable  to  high  energies. 
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IV.  PEL  Operating  Regimes 

there  are  three  regimes  of  free  electron  laser  operation.  The  operating 
regimes  can  be  characterized,  in  part,  by  the  velocity  distribution  relative 
to  the  phase  velocity  of  the  pondermotive  wave.  A  phase  space  plot  is  shown 
in  Figure  1. 

In  the  Raman  regime  the  beam  appears  cold  to  the  pondermotive  wave.  There 
are  no  particles  in  resonance  with  the  wave.  The  amplitude  of  the  wave  is  too 
small  to  trap  any  beam  particles  at  z«0.  As  the  wave  grows  in  space,  the 
amplitude  becomes  large  enough  to  trap  the  beam  particles,  thus  terminating 
the  linear  growth  phase. 

In  the  kinetic  Compton  regime,  the  beam  appears  warm  to  the  pondermotive 
wave.  There  are  particles  in  resonance  with  the  pondermotive  wave.  The 
linear  growth  is  proportional  to  the  slope  of  the  distribution  function.  This 
is  a  Landau  type  growth  mechanism.  Kinetic  theory  is  required  to  describe  the 
FEL  process  in  this  regime.  The  pondermotive  wave  is  resonant  with  a  range  of 
particles,  as  indicated  by  the  dashed  lines  in  the  figure.  The  growth 
mechanism  of  the  FEL  in  this  regime  has  a  quantvm  mechanical  analogy  with 
Compton  scattering.  It  is  not  a  single  particle  effect,  since  it  requires  a 
kinetic  description,  and  so  it  is  referred  to  as  kinetic  Conqpton.  The  phase 
space  plot  shows  the  growth  of  the  wave  in  space,  initially  the  wave 
aoplitude  is  small  compared  to  the  width  of  the  resonance.  The  pondermotive 
wave  will  grow  until  the  amplitude  becomes  large  enough  to  trap  the  resonant 
particles,  thus  terminating  the  linear  growth  phase. 


Raman 


Kinetic 

Compton 


Compton 


Figure  1 


VL\  -\  A 


Velocity  Distribution 


Phase  Space 


Velocity  distributions  and  phase  space  plots  of  the  three  FEL 


operating  regimes. 


In  the  Oompton  regime,  the  beam  energy  is  high  and  the  current  low. 

Single  particle  effects  dominate  over  collective  effects.  The  phase  velocity 
of  the  ponderaotive  wave  is  less  than  the  beam  velocity.  There  are  two  modes 
in  which  a  single  particle  Compton  FEL  may  be  operated.  If  the  amplitude  of 
the  ponderaotive  wave  is  too  small  initially  to  trap  the  beam  particles,  as  in 
the  Stanford  experiment,  the  gain  is  due  to  the  perturbation  of  beam  particles 
by  the  ponderaotive  wave.  The  interchange  between  wave  and  particle  energy  is 
oscillatory  at  the  same  period  as  the  ponderaotive  wave.  These  particles  are 
represented  by  a  shaded  region  above  the  ponderaotive  wave  ellipses  in  Fig.  1. 


If  the  wave  amplitude  is  sufficiently  large,  as  in  the  Los  Alamos 
experiment  (Warren  et.  al. ,  1982) ,  some  or  all  of  the  beam  particles  may  be 
trapped  in  the  potential  well  of  the  ponderaotive  wave.  The  gain  of  the  FEL 
in  this  region  is  due  to  the  loss  of  particle  energy  as  the  particles  rotate 
in  the  potential  troughs  of  the  ponderaotive  wave.  This  occurs  at  the  trapped 
particle  bounce  frequency.  The  maximum  gain  is  obtained  when  the  beam 
particles  are  at  the  bottom  of  the  ellipse  in  phase  space. 


The  efficiency  of  the  free  electron  laser  can  be  increased  in  all  of  these 
schemes  by  decreasing  the  phase  velocity  of  the  ponderaotive  wave  after  it  has 
trapped  the  beam  particles.  The  decrease  in  particle  energy  appears  as  an 
increase  in  wave  energy. 


I 


To  operate  a  high  current  beam  in  the  Raman  regime  requires 
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At  a  currant  density  of  1  kA/cm2  and  S  -  0.3  the  wavelength  must  be 
greater  than  0.3  cm  to  operate  in  the  Raman  regime  when  y  ®  3.  To  operate  a 
1  micron  PEL  in  the  Raman  regime  for  such  a  beam  vrould  require  y  ■  2600,  or 
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1.3  GeV,  and  a  wiggler  with  a  wavelength  of  135  meters.  Consequently ,  the 
Raman  FEL  is  not  viable  in  the  IR  region  unless  S  can  be  greatly  reduced. 


It  is  clear  from  Eq.  14  that  there  are  three  parameters  one  can  improve 
upon  to  extend  the  Raman  region.  The  kinetic  energy  of  the  beam  can  be 
increased/  the  current  density  decreased  or  the  scale  factor  in  the  emittance 
relation  decreased.  However,  there  is  a  limit  to  the  extent  one  can  increase 
the  energy,  decrease  the  current  density  and  maintain  a  collective 
interaction.  When  the  number  of  particle  in  the  debye  sphere  is  small,  the 
system  will  not  support  collective  oscillations. 


The  wavelength  condition  of  Bq.  14  is  most  sensitive  to  S,  so  it  is 
inportant  to  analyze  the  degree  to  which  S  can  be  reduced.  The  Lawson- Fenner 
relation  is  not  derived  from  first  principles  but  is  a  phenomenological 
relation  which  has  been  shown  to  hold  in  many  accelerators.  The  success  of 
the  Lawson-Fenner  relation  is  related  to  the  fact  that  the  reliable  current 
density  from  oxide  thermionic  cathodes  is  generally  <  10  A/cm  .  With  higher 
current  density  cathodes,  the  scale  factor,  S,  can  be  reduced  considerably. 

In  a  recent  experiment  using  cold  graphite  cathodes,  an  S  parameter  of  0.12 
was  obtained  with  a  14  kA  beam  by  aperturing  the  beam  to  about  30  percent  of 
the  diode  current  (Sloan  et  al,  1982) .  An  S  parameter  of  .13  at  4  kA  was 
obtained  by  aperturing  to  about  25%  of  the  current  from  a  spiral  knife-edge 
cold  cathode  in  a  20  MeV  induction  linac  (Kulke  et  al. ,  1981).  Finally,  a 
much  reduced  value  of  S  «*  0.04  at  1.5  kA  was  calculated  for  a  recent 
experiment  in  which  only  ~  5-10%  of  the  current  from  a  graphite  cold  cathode 
was  extracted  from  the  diode  (Parker  et  al.,  1982). 


The  S  parameter  is  best  understood  in  terms  of  the  brightness,  another 


invariant  that  is  related  to  the  emittance? 


B 

n 


ir2e  2 
n 


1  J 

2  2  "  * 
iTS  2*AY/Y 


(15) 


Thus,  the  brightness  is  a  measure  of  our  ability  to  make  a  high  current 
density,  cold  beam.  To  the  extent  that  the  emittance  is  an  invariant,  the 
brightness  is  determined  at  the  diode;  hence  if  we  congress  the  beam  at  a 
later  stage  we  increase  Ay/Y*  Therefore,  for  high  power  operation  it  is 
desirable  to  have  a  high  current  density  cathode.  From  the  above  three 
experiments  we  can  conclude  that  the  brightness  of  cold  cathodes  with 
appropriate  aperturing  is  a  factor  of  six  to  sixty  higher  than  that  of  hot 
cathodes.  This  will  extend  the  Raman  wavelength  region  by  a  corresponding 
amount.  However,  the  price  one  pays  for  cold  cathode  operation  is  diode 
closure.  This  can  limit  the  pulse  length  to  a  few  hundred  nanoseconds.  In 
addition,  the  efficiency  of  the  system  is  reduced  by  a  factor  of  5  to  10  when 
a  large  fraction  of  the  beam  is  terminated  by  aperturing. 

The  emittance  is  often  the  dominant  contributor  to  energy  spread. 

However,  the  self  potential  and  the  wiggler  gradient  can  also  contribute  to 
the  beam  energy  spread,  which  in  general  is  given  by 


where  the  second  term  is  the  contribution  from  the  self  potential  and  the 
third  is  from  the  wiggler  gradient*  From  this  we  can  find  an  optimum  radius 
of 
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For  high  current  beams  propagated  in  a  strong  axial  magnetic  field,  the 
equilibrium  radius  (Neil,  1977)  should  be  compatible  with  the  optimum  radius. 

The  variation  in  y  as  determined  by  the  emittance,  wiggler  gradients  and 
self  field  of  the  beam  results  in  a  velocity  spread  of  the  beam  that  is 
independent  of  time.  If  the  diode  voltage  is  time  dependent  (as  for  example 
due  to  diode  closure) ,  this  will  shift  the  output  wavelength.  If  the  gain  is 
sufficiently  high  in  an  oscillator  the  time  variation  will  cause  a  broad  band 
output.  The  time  variation  in  an  anplifier  experiment  can  cause  a  loss  of 
resonance  with  the  input  signal.  The  range  of  unstable  wave  numbers  is 
(Sprangle  et.  al. ,  1979) 
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where  k^  is  the  spatial  growth  rate.  TCiis  can  be  related  to  the  variation  in 
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Y  through  X  ■  X  /2y  .  Then  — I  <  4  k.  A  is  required.  Substituting 


the  maximum  growth  rate  value  for  k^ 
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Therefore,  the  requirements  on  the  time  variation  of  the  beam  energy  to 

keep  the  growth  of  the  output  wavelength  in  the  unstable  region  is 

3  3 

proportional  to  X  .  For  a  1  kA/cm  ,  1  MeV  beam  in  a  1  kG  wiggler  field  an 
output  wavelength  of  0.7  mm  requires  a  Ay/y  time  variation  of  less  than  1%. 

V.  Accelerator  Development 

During  the  past  fifty  years  accelerators  have  developed  primarily  along 
two  directions:  (1)  low  current,  high  voltage  and  (2)  high  current,  low 
voltage  devices.  Consequently  FEL  research  has  followed  similar  lines.  The 
history  of  traditional  accelerators  is  quite  well  known.  The  average  current 
in  these  accelerators  is  generally  less  than  an  ampere  with  micro  bunch 
currents  reaching  tens  of  amperes.  The  energies  acheived  in  electron 
linacs  exceeds  10 10  eV  and  the  beam  energy  in  a  proton  synchrotron  is 
approaching  10 12  eV.  High  current  accelerators  have  been  developed  primarily 
as  Marx  generators  driving  a  pulse  forming  line,  which  is  connected  to  a 
diode.  The  energy  of  these  devices  is  typically  around  1  MeV,  although  10  MeV 
devices  have  been  built.  The  pulse  current  from  a  pulse-power  line  diode  has 
reached  10  MA.  These  devices  do  not  scale  well  in  energy  because  the  total 
energy  is  achieved  in  a  single  gap,  hence  the  stress  on  the  insulator  becomes 
excessive  at  high  voltages.  A  free  electron  laser  requires  high  voltage  for 
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short  wavelength  operation  and  high  current  £or  efficiency  and  high  gain. 

Only  the  induction  linac  has  been  a  serious  candidate  as  a  FEL  driver  in  this 
parameter  range.  Experiments  at  Lawrence  Livermore  Laboratory  on  the 
Experimental  Test  Accelerator  (ETA)  have  acheived  currents  approaching  10  kA 
and  beam  energies  of  5  MeV.  A  scaled  up  version  of  this  accelerator  is 
designed  for  50  MeV,  10  kA,  30  ns  operation,  ftie  current  from  this  device  is 
probably  excessive  for  efficient  FEL  operation  and  the  pulse  length  too  short. 

A  long  pulse  induction  linac  (T  *  2  usee)  was  built  at  the  National 
Bureau  of  Standards  (Leiss,  et  al.,  1980)  and  is  currently  in  operation  at  the 
Naval  Research  Laboratory.  Figure  2  is  a  schematic  of  the  Long  Pulse 
Induction  Linac.  ttie  accelerator  consists  of  two  major  components;  (1)  an 
injector  and  (2)  an  induction  accelerator  module. 

The  electron  gun  in  the  injector  has  a  16.5  cm  diameter  tungsten  dispenser 
thermionic  cathode.  Electrons  are  accelerated  in  the  gun  through  a  series  of 
12  annular  electrodes,  spaced  by  ceramic  insulator  rings.  The  last  electrode 
supports  a  95%  transmission  tungsten  mesh  at  ground  potential.  The  gun  has 
also  been  operated  with  a  graphite  brush  cold  cathode  (Ramirez  and  Cook,  1980; 
Prohaska  and  Fisher,  1982) . 

The  electron  gun  is  immersed  in  an  oil  filled  tank.  The  gun  voltage  is 
fed  from  a  pulse  line  driving  a  12:1  step  up  transformer.  The  injector 
typically  produces  a  0.8  kA  beam  pulse  of  400  keV.  The  electron  beam  is 
transported  to  the  induction  accelerator  module  by  a  series  of  focusing  coils. 


INDUCTION  LINAC 


Figure  2.  Schematic  of  Long  Pulse  Induction  Linear  Accelerator. 


The  induction  accelerator  module  consists  of  two  core  sets.  One  of  the 
core  sets  give  a  4  to  1  step  up  voltage  and  the  other  a  5  to  1.  The  cores  are 
wound  with  0.001  inch  mild  steel  foil,  separated  by  0.00025  inch  mylar 
sheets.  Each  accelerating  gap  produces  about  200  kV  of  acceleration. 

Typically,  the  output  energy  of  the  electron  beam  generated  by  the  linac 
is  0.8  MeV,  the  current  approximately  0.8  kA  and  the  pulse  length  2  nsec. 
However,  for  the  data  reported  in  this  work  one  of  the  accelerating  gaps  was 
not  operating,  so  the  beam  voltage  was  0.55  MeV.  Also,  a  graphite  brush  cold 
cathode  was  used  and  the  0.8  kA  diode  current  was  reduced  to  0.2  kA  due  to 
losses  in  beam  transport.  The  temporal  variation  in  the  voltage  with  the  hot 
cathode  is  less  than  3  percent  over  1.6  usee.  This  is  the  longest  pulse 
induction  linac  in  existance,  with  a  pulse  duration  more  than  an  order  of 
magntiude  longer  than  other  induction  linacs.  Hie  pulse  length  becomes  an 
important  consideration  for  free  electron  laser  experiemnts  where  one  wishes 
to  study  the  nonlinear  dynamics  of  the  beam  or  efficiency  enhancement 
schemes.  Applications  that  require  a  significant  amount  of  energy  in  the 
radiation  field  also  require  long  pulses  to  avoid  the  problems  encountered 
with  excessive  electric  field  strengths  at  short  pulses. 

Figure  3a  shows  the  voltage  and  current  traces  of  the  injector  for  the 
long  pulse  induction  linac.  The  voltage  remains  constant  to  within  a  few 
percent  over  the  2  microsecond  pulse  length.  In  contrast  Figure  3b  shows  the 
voltage  and  current  from  the  graphite  cathode  diode  of  the  Febetron,  a 
relatively  long  pulse  Marx  type  generator.  The  oscillations  on  the  voltage 
and  current  are  due  to  the  finite  element  pulse  line  of  the  Febetron  Marx. 

The  diode  voltage  decreases  25  percent  in  250  nsec  as  a  result  of  diode 
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closure.  This  closure-induced  voltage  collapse  is  typical  of  cold,  plasma 


cathodes,  and  in  light  of  the  condition  of  Eg.  19,  is  a  key  problem  for  long 
pulse  FEL  operation.  Efforts  to  prevent  diode  closure,  such  as  magnetic 
insulation,  generally  result  in  higher  beam  emittance. 

The  relationship  between  the  cathode  plasma  velocity  and  the  diode 
parameters  can  be  determined  by  considering  the  simple  circuit  shown  in  Fig. 
4.  There  is  a  voltage  source  Vc,  an  internal  line  impedance  Z^,  and  a  time 
dependent  diode  impedance  zQ(t) .  The  current  is 


I  -  VC/(ZL  +  ZD>. 


(20) 


The  diode  ''It age  is 


VD  -  IZD  *  VcV(ZL  +  ZD> 


(21) 


and  the  diode  impedance  is 
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where  a  is  a  constant  and  d  the  anode  cathode  spacing,  we  will  assume 


d-d  -  (St  (23) 

o 

where  a  is  the  plasma  closure  velocity.  The  closure  velocity  for  high  current 
-told  cathodes  is  typically  ~  2.5  cm/ys.  Now  if  we  differentiate  Eq.  (21) 


Idealized  circuit  diagram  of  an  electron  beam  generator.  Source 
voltage  Vc  appears  across  line  impedance  Z L  in  series  with  diode 
impedance  Zn. 
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Hence  for  a  diode  that  is  initially  matched  we  can  obtain  a  simple  expression 
for  the  time  dependence  as  a  function  of  diode  spacing. 

Figure  5  is  a  plot  of  the  diode  voltage  as  a  function  of  time  for  Z^  -  ZD, 
ZL  *  10  ZD  and  zL  *  0.1  Zj>  when  t  ■  o»  When  the  line  impedance  is  ten  times 
the  diode  impedance,  the  device  tries  to  behave  as  a  constant  current  source 
and  the  voltage  decays  faster  than  in  the  matched  load  case.  When  the  line 
inpedance  is  small  compared  to  the  diode  impedance,  the  device  tries  to  act  as 


a  constant  voltage  source  at  the  expense  of  an  increased  current.  However 
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Figuxe  5.  Calculated  diode  voltage  and  current  for  three  different  diode 

impedances.  The  line  impedance  is  and  ZQ  is  the  diode  impedance 
at  t  *  0. 
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for  high  current  beams,  the  increasing  self  potential  of  the  been  with  current 
will  cause  the  beam's  kinetic  energy  to  decrease  faster  than  the  diode 
voltage. 

Although  diode  closure  in  high  current  cold,  cathode  machines  is  a 
limiting  factor  for  long  pulse  operation,  there  are  many  operational 
advantages  to  cold  cathodes.  Thus,  we  have  temporarily  replaced  the  hot 
cathode  in  the  induction  linac  with  a  graphite  brush  cold  cathode  of  the  sane 
area.  The  average  electric  field  at  the  cathode  surface  is  about  30  kV/cm. 
the  cathode  rapit’W  turns  on,  and  the  closure  problem  does  not  sens  to  be  as 
severe  as  in  the  Pebetron.  Perhaps  the  large  area,  low  current  density 
reduces  the  closure  velocity,  the  amittance  is  worse  than  with  the  hot 
cathode,  and  more  of  the  beam  is  lost  in  transport.  However,  the  emittance  of 
the  beam  transported  to  the  wiggler  appears  to  be  about  the  same  as  with  the 
hot  cathode. 

VI.  the  PEL  Apparatus 

Figure  6  is  a  sketch  of  the  free  electron  laser  experiment,  there  is  a 
uniform  axial  field  of  120  cm  length  which  is  varied  from  1  to  5  kG,  with  a 
typical  operating  field  of  2  kG.  the  beam  from  the  induction  linac  is  focused 
into  the  solenoid,  the  experiment  has  not  yet  been  run  with  the  mirrors  in 
place  as  shown,  but  operated  in  a  superradiant  mode. 

Two  wiggler  configurations  have  been  investigated:  (1)  a  pulsed  linear 
wiggler  and  (2)  a  radially  symmetric  diffusive  wiggler.  the  amplitude  of  the 
linear  wiggler  field  can  be  varied  from  0.1  kG  to  1  kG.  the  wiggler 


wavelength  in  3.0  cm,  and  the  overall  length  of  the  wiggler  is  120  cm.  The 
wiggler  amplitude  rises  adiabatically  in  30  cm,  has  a  uniform  straight  section 
of  60  cm  and  decays  adiabatically  in  the  last  30  cm. 

Two  diffusive  wigglers  have  been  used:  a  6  cm  period  aluminum  wiggler  and 
a  4.5  cm  copper  wiggler.  The  perturbed  axial  component  of  the  on-axis  field 
is  4  percent  of  the  axial  field  for  the  6  cm  wiggler  and  6  percent  for  the  4.5 
cm  wiggler. 

The  output  radiation  is  extracted  through  a  large  area  window,  and  the 
power  and  spectrum  are  determined  by  gas  breakdown  thresholds  or  with  high 
pass  filters  and  calibrated  detectors. 

A.  Linear  Wigglers 

There  are  several  advantages  to  linear  wigglers,  including  ease  of 
assembly,  changing  the  periodicity  to  operate  at  different  wavelengths,  and 
tapering  the  period  and/or  field  anplitude  for  efficiency  enhancement 
(Sprangle  et  al. ,  1979).  As  a  result  of  our  experiments  and  analysis,  we  find 
that  a  linear  wiggler  in  an  axial  guide  field  produces  an  elliptical 
polarization.  The  ratio  of  the  major  to  minor  axis  of  the  ellipse  is  kwv2/^0/ 
where  kw  is  the  wiggler  wave  nvanber,  vz  the  axial  beam  velocity,  and  the 
cyclotron  frequency  in  the  axial  guide  field.  As  a  result  of  the  asymmetry  of 
the  wiggler,  there  are  no  focusing  forces  in  one  of  the  directions 
perpendicular  to  the  beam  propagation  and  an  electron  drift  results.  Although 
this  drift  can  be  very  small  in  some  parameter  regimes,  it  is  always  present 
in  linear  wigglers  immersed  in  an  axial  guide  field  unless  additional  focusing 


forces  are  provided.  This  problem  was  avoided  in  a  high  current  device  using 
a  linear  wiggler  by  eliminating  the  solenoidal  field  (Phillips,  1960). 

Focusing  was  provided  by  increasing  the  wiggler  field  at  the  edges. 

To  understand  how  the  elliptical  polarization  and  perpendicular  drift 
arise,  we  will  solve  the  equations  of  motion  in  the  combined  wiggler  and  axial 
fields.  The  linear  wiggler  field  components  are 


B  »  b(z)  B  cosh  k  y  cos  k  z  (28) 

wy  www 


B  =■  -  b(z)  B  sinh  k  y  sin  k  z 
wz  w  w^  w 

where  Bw  is  the  peak  wiggler  field  on  axis  and  b(z)  is  an  adiabatic  taper  in 
the  wiggler  anplitude  at  the  entrance  and  exit  of  the  wiggler.  The  amplitude 
increases  to  its  full  value  in  ten  periods  in  our  experiment. 


We  define  =  qBQ/Ya  and  *  qBw/Ym.  To  simplify  the  analysis  we  will 
make  the  following  assumptions 


B  v  v 

ky  <  <  1 
B  v  v  vr 

o  z  z 


and  v  *  v.  “  constant.  The  single  particle  equation  of  motion  then  gives 

Z  D 
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v  *  |S1  v  -ft  cosh  k  y  cos  k  z  v  ) 
x  o  y  w  w  w  z 


v  -  -v  ft 
y  x  o 


v  -  v  ft  cosh  k  y  cos  k  z 
z  x  w  w  w 


where  we  have  neglected  the  correction  to  B2  from  the  wiggler  since  B^Bjj  <<  1 


and  v  <<  v  . 
y  z 


We  take  the  derivative  of  each  of  the  components  with  respect  to  time  and 


use  the  assumptions  of  Eg.  29  to  get 


“2  2 
v  +  ft  v  -  ft  k  cosh  k  y  sin  k  z  v 
xoxww  w  wz 


Similarly 


v  +  ft  v  ■  ft  ft  cosh  k  y  cos  k  z  v 
yoyow  w  wz 


We  have  the  following  solutions 


ft  k  cosh  k  y  sin  k  z  v 
w _ jr _ w  z 

.  2  ,2  2 
ft  -  k  v 
o  w  z 


ft  ft  cosh  k  v  cos  k  z  v 
o  w _ w_ _ w  2 

n  2  ,2  2 

ft  -  k  v 
o  w  z 


ft  v  cosh  k  y  cos  k  z 
w  z  w^  w 
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ft  ft  cosh  k  y  sin  k  z 
ow _ w _ w 


k  (ft  ' 
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The  displacement,  in  the  x  direction  is  proportional  to  the  wiggler  field 
and  axial  velocity,  The  displacement  in  the  y  direction  occurs  only  if  there 
is  a  nonzero  solenoidal  magnetic  field.  It  is  proportional  to  both  the 
solenoidal  and  Higgler  field.  This  displacement  is  a  result  of  the  Lorentz 
force  arising  from  the  x-component  of  vw  and  BQ . 


The  ratio  of  the  maximum  excursion  in  the  x-direction  to  that  in  the  y- 
direction  is  seen  from  Eq.  33  to  be 


Ax/ Ay  ■  k  v  /ft  .  (34) 

w  z  o 

Thus,  this  ratio  is  a  measure  of  how  close  ft  is  to  the  cyclotron  resonance 

o 

which  occurs  at  ft  ■  k  v  .  Since  k_v_  is  approximately  constant  in  the 

O  w  z  w  z 

experiment,  the  ratio  of  the  maximum  displacements  Ax/Ay  should  vary 

as  1/ft  . 
o 

Figure  7  is  a  series  of  exposures  made  from  x-rays  produced  when  the  beam 
strikes  a  target  placed  beyond  the  wiggler.  These  photos  give  us  the  shape  of 
the  beam  versus  solenoidal  field.  In  Fig.  7(a),  the  wiggler  field  is  zero  and 
the  beam  is  approximately  circular.  In  Fig.  7(b),  with  a  500  Gauss  wiggler 
field  the  beam  becomes  elongated  in  the  x-direction.  The  linear  wiggler  field 
is  in  the  y  direction,  consistent  with  the  above  analysis.  As  the  solenoidal 
field  is  increased  with  the  wiggler  field  constant,  the  shape  goes  from 
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B0  -  2.6  Byy  -  0.5  B0  =  3.48  Bw  -  0.5 

(0  (d) 

Figure  7.  Time  integrated  x-ray  exposure  of  beam  striking  a  target  placed 


beyond  linear  wiggler  (wiggler  field  is  vertical)  for  various 
values  of  solenoidal  field. 
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elliptical  to  circular,  as  shown  in  Figs.  7(b)  -  7(d).  Hence,  we  would  expect 
a  nearly  linearly  polarized  output  for  (b)  and  a  nearly  circular  polarization 
for  (d) . 

Figure  8  is  a  plot  of  the  ratio  Ax/A y  as  a  function  of  the  solenoidal 
field.  The  results  are  plotted  three  ways.  The  solid  line  is  the  theory,  the 
small  dots  were  confuted  numerically  and  the  circled  x  points  are 
experimental.  At  small  values  of  magnetic  field  the  approximations  break  down 
and  there  is  a  large  variation  in  the  results. 

We  have  demonstrated  that  the  polarization  of  a  linear  wiggler  in  a 
solenoidal  magnetic  field  may  be  varied  by  changing  the  magnitude  of  the 
solenoidal  field.  In  addition,  we  see  from  Eqs.  34  that  a  change  in  the 
amplitude  of  the  wiggler  field  alone  does  not  change  the  polarization.  We 
have  confirmed  this  by  propagating  a  beam  the  entire  length  of  the  wiggler  in 
a  solenoidal  magnetic  field  of  2  kG.  The  wiggler  field  is  varied  up  to  1  kG 
without  changing  the  shape  or  intensity  of  the  beam. 

As  a  result  of  the  asymmetry  of  the  wiggler,  there  are  no  focusing  forces 

in  the  x-direction.  This  results  in  a  net  drift  of  off  axis  particles  (Pasour 

et  al.,  1982).  This  drift  is  shown  clearly  in  Fig.  9,  which  consists  of 

electon  trajectories  in  the  x-y  plane.  The  trajectories  are  calculated  by 

numerically  solving  the  equations  of  motion,  including  self-fields.  Initially 
■¥ 

the  electrons  E  x  B  drift  azimuthally  as  they  are  injected  into  to  the 
adiabatically  increasing  wiggler  field.  However,  when  they  reach  the  constant 
amplitude  portion  of  the  wiggler,  the  electrons  drift  rapidly  in  a  direction 
perpendicular  to  both  the  axial  and  wiggler  fields. 
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AX /AY 


•  Computed  Numerically 
0  Experimental 


B(kG) 

Figure  8.  Plot  of  beam  elongation  (dx/Ay)  vs.  axial  magnetic  field  in  a 
linear  wiggler  (Bw  in  y-direction) .  Solid  line  is  theoretical 
curve  (Eq.  34) ,  dots  are  values  obtained  by  numerical  integration 
of  orbit  equations,  and  crosses  are  experimental  results. 


Figure  9.  Electron  trajectories  in  a  linear  wiggler  and  an  axial  magnetic 


field.  In  this  case,  B  =  1  kG,  B  =  250  A,  X  =3  cm,  and 

W  0  w 

Y  =  2.2. 
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An  approximate,  empirical  expression  for  the  drift  has  been  found  which  is 
in  good  agreement  with  results  from  coirputer  calculations  of  the  electron 
orbits  when  nw  <  0o  <  kv2*  In  terms  of  the  wiggler  gradient,  the  expression 
is 


2  ft 


o 


V(B  >2  B 
w  o 

~~r~2 x  r  ' 

2(  B  )  o 

w 


(35) 


2  2  2 

where  v.  =  v  +  v  .  This  expression  has  the  same  form  as  the  usual  grad  B 
I  x  y 

drift  and  can  be  thought  of  as  arising  from  the  crossed  wiggler  gradient  and 
axial  field.  It  can  also  be  written  as 


v.  = 


2  vz  ft 


k  v  (ft  k  v  )2  cosh  k  y  sinh  k  y 
w  z  vwwz'  vr  w 


(36) 


(ft  -  k  v  2) 
^  o  w  z  ’ 


This  expression  is  reasonably  valid  for  kwy  <_  0.8,  Bwy  >>  Bwz,  2uid 

<<  v^.  Note  that  the  drift  is  very  small  near  the  axis  but  increases 
exponentially  with  kwy. 

Table  1  compares  drift  velocities  for  various  cases  as  calculated  by  the 
computer  trajectory  code  to  those  calculated  from  the  Bg.  36.  In  general  the 
agreement  is  very  good.  One  consequence  of  the  large  drift  near  resonance  is 
that  it  limits  the  degree  of  gain  enhancement  achievable  through  the  magneto¬ 


resonance  effect 


1  ABLE  I 


Comparison  of  drift  velocities  computed  from  the  code  with  those  from  Eq.  9  for 
various  cases. 
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2.2 

0.29 

0.4 

0.047 

0  .050 
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2.2 

0.57 

0.4 

0.042 

0  .046 
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0.5 

2.2 

0.29 

0.4 

0.011 

0  .013 
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1 

2.2 

0.29 

0.2 

0.019 

0  .022 

2 

1 

3.0 

0.21 

0.4 

0.029 

0.034 

4  - 

1 

3.0 

0.42 

0.4 

0.019 

0  .022 
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0.5 

3.0 

0.42 

0.4 

0.0051 

0  .0056 

4 

1 

10.0 

0.11 

0.4 

0.0043 

0.0049 

10 

5 

10.0 

0.28 

0.4 

0.057 

0.056 
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B.  Helical  Higgler 


Because  the  linear  wiggler  drift  is  so  sensitive  to  slight  off-centering 
or  large  diameter  of  the  beam,  a  helical  wiggler  vrould  seem  to  be  somewhat 
better  suited  to  the  experiment,  at  least  initially  when  beam  focusing  is  not 
optimized.  If  the  trajectory  calculation  is  repeated  with  the  same  parameters 
as  those  in  Fig.  9  but  with  a  helical  wiggler  substituted  for  the  linear  one, 
it  is  found  that  the  electron  drifts  azimuthally  and  is  well  confined,  as 
shown  in  Fig.  10  (Pasour  et  al. ,  1982).  This  behavior  results  from  the 
radially  increasing  wiggler  field,  which  for  k^r  <0.8  can  be  written  as 
(Blewett  and  Chasman,  1977) 


B  -  B  (1  +  |  k  2r2)  sin  (0  -  k  z) 
r  w  8  w  w 

I 

B0  -  B  (1  +  ^  k  2r2)  cos  (9  -  k  z)  (37). 

0  w  8  w  w 


B  »  -  krB  (L  +  ~  k  2r2)  cos  (0  -  k  z). 
z  w  8  w  w 


Although  a  helical  wiggler  is  in  general  more  difficult  to  construct  than 
a  linear  one  and  is  not  nearly  as  amenable  to  tapering,  the  orbital  stability 
that  it  provides  is  a  major  advantage  when  an  axial  guide  field  is  used. 

Also,  it  is  possible  to  taper  the  amplitude  of  the  wiggler  field  in  the 
helical  wiggler  by  carefully  varying  the  winding  radius  as  a  function  of  z. 


C.  Diffusive  Wiggler 


An  axially  symmetric  wiggler  can  be  constructed  by  placing  conducting 
rings  in  a  pulsed  solenoidal  magnetic  field  (Jacobs  et  al. ,  1980).  As  the 
solenoidal  field  difffuses  in,  it  induces  eddy  currents  in  the  conducting 
rings  which  generate  an  opposing  field.  This  results  in  a  magnetic  field  with 
a  periodic  radial  conponent.  This  process  is  illustrated  in  the  computer 
plots  of  Fig.  11. 

The  magnetic  field  at  a  distance  r  from  the  the  axis  is  approximately 


♦ 

B  a  Z 


[B  +  B,I  (k  r)  sin  (k  z) ]  -  rB.I,  (k  r)  cos  (k  z)  , 
o  lo  w  w  1  1  w  w 


(38) 


where  IQ  and  1^  are  modified  Bessel  functions  and  B^  is  the  amplitude  of  the 
axial  field  modulation  on  axis. 

Figure  12  is  a  computer  plot  of  the  axial  conponent  of  the  magnetic  field 
of  the  diffusive  wiggler  used  in  the  system.  When  the  wiggler  is  not  present, 
the  solenoidal  field  is  2  kG.  Hence,  the  wiggler  not  only  modulates  the 
solenoidal  field,  it  reduces  the  average  value.  To  make  a  smooth  transition 
into  the  wiggler  field,  the  conducting  rings  extend  all  the  way  to  the  end  of 
the  solenoid  where  the  field  drops  to  one  half  the  peak  value. 

Since  the  diffusive  wiggler  is  axially  symmetric,  it  does  not  produce  the 
kind  of  radial  particle  drift  which  is  characteristic  of  linear  wigglers  in  an 
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Axial  Position  (cm) 

Figure  12.  Plot  of  B  (r=0)  vs.  z  in  the  diffusive  wiggler.  Face  of  solenoid 

2 


is  at  z  *  0. 


axial  field.  However,  the  perpendicular  component  of  the  wiggler  vanishes 
along  the  axis  of  the  diffusive  wiggler,  so  that  only  electrons  that  are 
significantly  off  axis  can  participate  in  the  type  of  FEL  interaction 
described  previously.  It  has  been  pointed  out,  however,  that  an  interaction 
that  relies  on  the  transverse  energy  in  the  beam  can  be  exploited  to  obtain 
radiation  from  such  a  wiggler.  This  device  has  been  called  the  Lowbitron 
(McMullin  and  Bekefi,  1981).  The  radiation  frequency  is  approximately 

u>  -  (1  +  0  2)  Y  2  (v  k  /y) .  (39) 

z  z  z  w  o 

This  radiation  is  a  result  of  electrons  interacting  with  the  axial 
component  of  the  wiggler  field.  'Bie  Lowbitron  interaction  requires 

<  1.  The  electron  gyroradius  must  be  sufficiently  small  that  the 
transverse  field  modulation  felt  by  the  electron  can  be  neglected  compared 
with  the  longitudinal  modulation. 

VII.  Radiation  Measurements 

We  have  carried  out  a  series  of  measurements  on  the  radiation  using  a 
diffusive  wiggler.  Radiation  measurements  have  primarily  consisted  of  a) 
spectral  analysis  using  cylindrical  cut-off  filters  and/or  a  gas  breakdown 
spectrometer  and  b)  approximate  power  measurements  using  either  pressure 
thresholds  for  gas  breakdown  or  calibrated  crystal  detectors  and 
attenuators.  Careful  consideration  has  been  given  to  the  identification  of 
the  interaction  modes,  both  cyclotron  and  FEL,  and  scaling  measurements  have 
been  performed  to  verify  these  modes. 
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A.  Gas  Breakdown  Spectrometer 


The  radiated  power  can  be  estimated  from  the  pressure  at  which  gas 
breakdown  occurs.  An  empirical  formula  for  the  rms  breakdown  electric  field 

®b  is 


Eb  «  AP[l  +  (a/PX)2l1/2,  (40) 

where  A  =  3000  VtoT^Torr”*  and  a  -  0.9  Torr-m  for  air  or  nitrogen,  P  is  the  gas 
pressure  in  'Dorr,  and  X  is  the  wavelength  in  meters.  Equation  (40)  has  been 
checked  against  McDonald's  gas  breakdown  data  (McDonald,  1966)  and  agrees  in 
the  worst  case  to  within  30%.  This  formula  is  valid  on  the  high  pressure  side 
of  the  breakdown  curve  and  for  X  <  30  cm,  pulse  lengths  >1  us,  low  repetition 
rate  (<  100  pps),  and  electron  diffusion  lengths  small  conpared  to  the  chamber 
dimensions. 


If  the  radiation  is  reflected  from  a  boundary,  a  standing  wave  is  set  up 
with  an  amplitude  given  by 

E2  -  Z*  (l  +  r2  -  2r  cos  [2  £  •  x])  (41) 


where  E^  is  the  magnitude  of  the  incident  wave,  T  the  reflection  coefficient, 

and  k  the  wavenumber.  The  distance  d  between  the  peaks  in  the  standing  wave 

♦  + 

is  determined  by  setting  2  k  •  x  -  2x,  so  that 


d 


X 

2  cos  9 


(42) 


where  X  is  the  wavelength  of  the  radiation  and  9  the  angle  of  incidence  with 
respect  to  the  surface  normal.  By  adjusting  the  gas  pressure  so  that  the 
standing  wave  field  amplitude  is  slightly  higher  than  the  breakdown  field, 
localized  gas  breakdown  will  occur  at  the  standing  wave  peaks.  Hien  X  can  be 
determined  by  measuring  the  distance  between  adjacent  breakdown  spots,  and  the 
power  can  be  estimated  from  the  pressure  required  to  initiate  observable 
breakdown. 

Figure  13  shows  a  schematic  of  the  spectrometer.  Radiation  is  collected 
in  the  horn  on  the  left.  It  then  travels  through  a  high  pass  filter  and 
expands  through  the  horn  on  the  right  to  a  collimating  or  weakly  focusing 
lens.  Hie  microwaves  are  reflected  from  a  metallic  boundary  located  inside 
the  gas  filled  chamber. 

Figure  14  is  a  time  integrated  photograph  of  the  resulting  interference 
pattern.  Microwave  radiation  has  entered  from  the  left  of  the  photograph 
where  the  lens  is  located  (but  cannot  be  seen)  and  is  reflected  off  a  copper 
plate  located  at  the  right  of  the  photograph  (outlined  by  a  light  ring).  All 
the  white  spots  are  due  to  light  produced  when  the  gas  breaks  down.  Type  57 
(ASA  3000)  Polaroid  film  was  used  with  a  Graphflex  camera  (f/4.5)  to  obtain 
these  results.  The  gas  density  was  selected  so  that  enough  light  was 
available  for  photographing  the  spots,  but  was  kept  low  enough  to  prevent 
microwave  reflections  from  the  plasma.  Thus  for  the  measured  pressure  (25 
Torr)  and  wavelength  (4.5  cm),  the  electric  field  is  determined  from  Bq.  40  to 
be  about  1  kV/cm.  From  the  spot  diameter  of  2.5  cm,  the  power  is  estimated 
for  a  plane  wave  to  be  8  kW. 
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This  technique  underestimates  the  power  because  there  is  considerable 


energy  outside  the  radius  of  the  observed  spot.  For  example,  in  the  case 
above  the  spot  diameter  is  nearly  a  factor  of  2  smaller  than  the  wavelength, 
so  clearly  the  radiation  envelope  is  larger  than  the  measured  spot  size. 

We  have  identified  this  radiation  as  a  TEQ1  mode  that  is  excited  at  the 
2nd  harmonic  of  the  electron  cyclotron  beam  interaction.  With  similar 
measurements  we  have  observed  more  than  100  lew  of  power  at  a  wavelength  of  7.5 
cm.  This  mode  was  identified  as  a  TE^  mode  excited  at  the  fundamental 
cyclotron  beam  interaction.  The  electron  beam  current  in  these  experiemnts 
was  130  A,  the  voltage  550  keV  for  2  usee,  and  the  wiggler  period  was  6  cm. 
The  gas  breakdown  spectrometer  has  been  described  in  more  detail  elsewhere 
(Mako  et  al.  ,  1982)  . 

B.  Mode  Analysis  and  Experimental  Be  suits 

The  strongest  interactions  are  expected  when  the  phase  velocities  of  the 
beam  modes  and  waveguide  modes  are  equal.  The  dispersion  relation  for  the 
waveguide  mode  is 


2  2.2  2  .... 

to  *  a)  +  k  c  (43) 

co 

where  u>  is  the  cutoff  frequency.  The  dispersion  relation  for  the  free 
co 

electron  laser  is 
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Hie  last  term  is  snail  at  the  current  densities  of  this  experiment  and  is 
neglected.  The  cyclotron  mode  dispersion  relation  is  given  by 

a)  =  k  v  +  n  Q  /y,  (45) 

20 

where  Q  =  eB^m  and  y  have  been  explicitly  written  to  make  the  energy 
scaling  clear. 

Then  the  intersections  of  the  FEL  and  cyclotron  modes  with  a  waveguide 
mode  are  found  from  Eqs.  43-45  to  be 
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where  w  *  2irc/aX  ,  a  is  the  waveguide  radius  and  X„„  =  3.41  for  TE-,  ,  ,  2.61 
co  nm  nm  11 

for  TMQ1,  2.06  for  TE2^,  1.14  for  TM^ ^  and  1.14  for  TMq2  modes. 

A  sketch  of  three  dispersion  relations  is  shown  in  Fig.  15.  The  cyclotron 

mode  will  be  above  the  FEL  mode  if  nQ  /y  >  k  v  .  There  are  typically  both 

o  w  z 

high  frequency  and  low  frequency  intersections  of  the  FEL  and  cyclotron  modes 
with  the  waveguide  mode.  However,  the  high  frequency  FEL  intersection  is  the 
one  which  results  in  the  usual  FEL  interaction.  Also,  the  gyrotron  typically 
operates  at  the  low  frequency  cyclotron  intersection,  but  a  cyclotron 
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interaction  is  also  possible  at  the  high  frequency  intersection.  Figure  16  is 

a  plot  of  the  waveguide  beam  mode  interactions  that  are  possible  for  the 

parameters  of  the  exper invent.  Both  the  first  and  second  harmonics  of  the 

electron  cyclotron  interaction  are  shown,  assuming  y  “  1.6.  The 

z 

justification  for  this  assumption  will  be  given  below.  The  arrows  in  the 
figure  indicate  the  frequencies  at  which  the  strongest  interactions  were 
observed  experimentally. 

To  positively  identify  the  FEL  mode  at  f  *  16  GHz,  we  have  varied  the 

magnetic  field  and  wiggler  wavelength.  It  turns  out  that  the  cyclotron  mode 

is  very  insensitive  to  BQ  in  our  parameter  range.  Because  the  magnetic  moment 

2 

is  approximately  conserved,  the  leading  term  for  the  cyclotron  wave  12  y  is 

o  z 

nearly  a  constant  for  small  changes  in  the  magnetic  field.  Thus,  the 

frequencies  of  the  electron  cyclotron  modes  increase  only  slightly  with 

increasing  magnetic  field.  However,  the  FEL  frequency  increases  with  a 

2 

decreasing  magnetic  field  because  Yz  varies  as  the  inverse  of  Bq. 

Figure  17  contains  the  results  of  the  mode  identification  experiments. 

The  measurements  were  made  using  high  pass  cutoff  filters,  so  each  bar  is  an 
indication  of  the  resolution.  Ibis  technique  integrates  the  total  power  above 
the  cut  off  frequency  of  the  particular  filter.  Thus,  the  power  in  a 
particular  band  is  just  the  power  measured  with  that  filter  minus  the  power 
measured  with  the  next  smaller  diameter  filter. 

The  arrow  is  the  calculated  FEL  frequency  with  y  *  1.6.  At  X  =  7  cm 

z  w 

and  Bq  -  2.45  kG  the  calculated  frequency  is  11.6  GHz  and  this  is  where  we 
see  a  factor  of  two  increase  in  the  signal  amplitude.  When  the  magnetic  field 
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Experimental 


1  1 


I 


I 


Frequency  (GHz) 


Figure  16.  Beam-waveguide  mode  interactions  which  are  possible  with  the 


experimental  parameters.  The  arrows  denote  the  measured  output 
frequencies . 


Frequency  (  GHz ) 


Figure  17.  Frequency  scaling  of  one  particular  mode  with  magnetic  field 
and  wiggler  period,  as  determined  by  using  a  series  of  high- 
pass  cut-off  filters.  Except  as  limited  by  the  detector  band¬ 
width,  the  signal  amplitude  shown  in  each  frequency  interval  is 
proportional  to  the  integrated  poweT  at  all  frequencies  above 
the  cut-off  of  the  corresponding  filter.  Thus,  power  in  a 
particular  interval  is  indicated  by  a  decrease  in  signal  amplitude 
at  the  next  higher  interval.  The  arrows  locate  the  theoretical 
frequency  of  the  TE^  FEL  mode,  which  agrees  well  with  the 
observed  spectra. 


V 


is  reduced  to  2.17  kG,  the  FEL  frequency  should  increase  to  13.2  GHz,  and  we 

observe  a  factor  of  two  increase  in  that  cutoff  filter.  At  X  -  6  cm  and  B  = 

w  o 

2.17  kG,  the  FEL  interaction  should  occur  at  16.5  GHz  and  we  observe  a  factor 

of  10  increase  in  that  cutoff  filter.  We  have  also  observed  that  the  electron 

cyclotron  mode  at  21  GHz  increases  with  magnetic  field.  We  have  looked  for 

2 

radiation  from  the  Lowbitron  interaction  at  oa  *  2y  (ft  /y  +  k  v  ) .  This 

1  z  o  1  w  z 

interaction  requires  perpendicular  velocity  and  interacts  with  variations  in 
the  z  cooponent  of  the  wiggler  field.  We  did  not  see  Lowbitron  radiation  that 
was  conparable  in  amplitude  with  the  FEL  and  2nd  harmonic  cyclotron  radiation. 

In  summary,  we  have  observed  free  electron  laser  radiaton  from  the  long 
pulse  induction  linac.  The  FEL  was  operated  in  the  superradiant  modo.  Hie 
power  radiated  in  the  FEL  mode  was  approximately  10  kW  with  the  beam 
interacting  with  the  TE^  waveguide  mode.  Second  harmonic  electron  cyclotron 
radiation  of  conparable  amplitude  was  observed  in  the  TEQ^  mode.  Most  of  the 
power  resulted  from  the  low  frequency  cyclotron  interaction.  More  than  100  kW 
of  radiated  power  at  the  fundamental  electron  cyclotron  mode  was  observed  in 
the  TE.^  mode.  In  all  cases  the  duration  of  the  radiation  was  2  usee.  The 
thermal  spread  of  the  beam  in  the  region  requires  FEL  operation  in  the  kinetic 
Coop  ton  regime. 

The  radiation  spectrum  is  very  sensitive  to  the  average  value  of  y  , 

z 

which  was  determined  by  fitting  the  observed  spectrum  with  interaction 
frequencies  calculated  using  an  assisted  value  of  <y^>.  Then  small  variations 
in  the  magnetic  and  wiggler  period  were  made  to  determine  if  the  mode  was  an 
FEL  or  cyclotron  mode.  In  this  manner  we  were  able  to  determine  that  <YZ> 


had  to  be  1.6  +  .2 


Another  way  of  estimating  y  is  by  measuring  the  beam  radius  in  the 

z 

uniform  field  of  the  solenoid.  Since  the  cathode  in  this  experiment  is  in  a 
magnetic  field  free  region,  the  canonical  angular  momentum  P  is  zero.  Then 


ymr  v0  *  erAQ, 


(48) 


where  Ag  is  the  magnetic  vector  potential  which  for  a  uniform  magnetic  field 

is  Aa  *  r  B  /2  .  Then 
a  z 


.  eB  r 
1  z 
2  ymc 


(49) 


Experimentally,  the  rms  beam  radius  was  measured  from  time  averaged  x-ray 

scintillation  pictures  to  be  1.3  cm  when  B_  *  2.17  kG  and  y  *  2.1.  Thus 

z 

conservation  of  PQ  alone  gives  B^“  0.4  at  this  rms  beam  radius.  Then 


Y/(l  +  Bj2  Y2) 


In  addition,  computer  simulations  of  the  induction  linac  free  electron 
laser  configuration  have  been  carried  out  (Thomson  et.  al. ,  1982).  The  mean 
value  of  Y  from  the  simulations  was  1.61.  Most  of  the  energy  spread  comes 
at  the  transition  between  the  induction  linac  transport  system  and  the 
solenoidal  field.  Figure  18  is  a  plot  of  the  induction  linac  FEZ,  magnetic 
field  profile.  The  diode  is  in  a  field  free  region  at  z  ■  0  and  the  last 
focusing  coil  of  the  linac  is  at  350  cm.  The  edge  of  the  solenoid  is  at 
approximately  375  cm.  The  2  kG  magnetic  field  is  required  for  a  beam 
equilibrium  radius  that  is  consistent  with  a  practical  wigglef  diameter. 


(kG) 
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Figure  19  is  a  plot  of  results  from  a  static  computer  simulation  of  the  beam 
transport  from  the  cathode  into  the  solenoid.  The  last  accelerating  gap  of 
the  induction  module  is  not  energized,  so  "y  =  2.1  as  the  beam  enters  the 
last  focusing  coil.  When  the  beam  has  reached  the  uniform  field  region  <y^> 
is  1.6. 

This  problem  of  beam  transport  and  matching  represents  one  of  the  major 

design  concerns  for  high  current  free  electron  lasers.  The  equilibrium  of 

high  current  beams  is  sensitive  to  self  electric  and  magnetic  fields  of  the 

beam.  Abrupt  changes  in  the  wall  diameter  or  magnetic  field  profile  induces 

betatron,  or  zero  frequency  cyclotron  oscillations  at  u>  *  0  *  kf?c  -  Q  /y ,  so 

o 

the  period  is  X  *  2ir8cy/fl  • 

o 

A  spread  in  B  across  the  beam  radius  will  cause  these  oscillations  to  phase 
mix  after  a  few  periods.  The  perturbation  then  appears  as  an  effective 
temperature. 

In  the  present  simulation  the  initial  normalized  emittance  was 

e  *  180  it  mrad-cm,  whereas  the  effective  final  emittance  was 
n 

e  *  yrB.  _  *  940  ir  mrad-cm.  Almost  all  of  this  increase  in  emittance 

n  i  rms 

came  in  the  transition  to  the  solenoid.  This  is  why  the  final  emittance  for 
the  hot  and  cold  cathode  is  about  the  same.  These  problems  can  and  have  been 
avoided  in  Marx-pulse  line  beam  generators  by  immersing  the  diode  in  a  very 
strong  magnetic  field.  This  minimizes  the  radial  excursion  and  effective 
emittance  growth.  This  approach  becomes  difficult  when  long  pulse  times  are 
required  and  one  must  use  large  area  hot  cathodes  to  get  kiloanqpere  electron 
beams.  One  can  conclude  from  this  that  the  design  of  the  beam  transport 
system  for  an  induction  linac  FEL  is  at  least  as  important  as  the  diode 
design. 
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VIII.  Conclusions  and  Future  Directions 

Improvements  in  the  beam  transport  are  clearly  needed  to  operate  in  the 
Raman  regime.  One  of  the  advantages  of  operating  in  the  cold  beam  limit  is 
that  the  efficiency  of  the  PEL  is  higher  and  can  be  enhanced  by  tapering  the 
wiggler  period.  Once  the  beam  is  trapped  in  the  pondermotive  wave,  the  phase 
velocity  of  the  wave  can  be  reduced  gradually  to  extract  energy  from  the 
beam.  Figure  20  is  a  plot  of  efficiency  vs  distance  for  an  untapered  and 
tapered  wiggler.  These  results  are  from  a  computer  simulation  of  the  present 
experiment.  Whereas  the  intrinsic  efficiency  is  about  5  percent,  the  tapered 
wiggler  efficiency  is  in  excess  of  20  percent. 


To  go  to  shorter  wavelengths  the  FEL  can  be  operated  as  a  two  stage 

2 

device.  The  output  from  stage  one  at  X  ■  Xw/2y  can  be  used  as  the  output 
for  the  second  stage  to  give 


X 


2 


(51) 


To  obtain  a  hi<gfri  voltage  high  current  accelerator  and  achieve  a  high  gain, 
high  power  short  wavelength  FEL,  the  long  pulse  induction  module  can  be 
converted  into  a  racetrack  accelerator  (Roberson,  1981?  MondelLi  and  Roberson, 
1982).  This  takes  advantage  of  the  long  pulse  to  give  a  voltage 


V  -  V  T/t  (52) 

9 

where  Vg  is  the  voltage  gain  of  the  module,  T  the  time  the  module  is  on,  and 
x  the  time  it  takes  the  beam  to  go  around  the  racetrack.  In  the  present 
case  a  30  nsec  path  length  would  result  in  a  26  MeV  beam  with  the  2  vs  ,  0.4 


MV  module 


0  30  60  90  120  150  100  210  240 

WIGGLER  LENGTH  (CM) 

Figure  20.  Theoretical  FEL  efficiency  vs.  interaction  length  with  an  without 
a  tapered  wiggler.  The  parameters  used  for  the  calculation  are 
those  of  the  experiment. 
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